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ABSTRACT 

Context. When modeling the density and grain size distribution in debris disks, the minimum particle size is often significantly larger 
than the corresponding blowout size. While the dust particles are usually modeled as compact, homogenous spheres, we instead in- 
vestigate the impact of porosity. 

Aims. The optical properties of porous particles are determined, and the influences of porosity on the blowout size and dust tempera- 
tures investigated. 

Mettiods. Using the method of discrete dipole approximation, we calculate the scattering and absorption cross sections of porous 
particles and derive the blowout size and the behavior of the dust temperature. 

Results. We determine the elficiency factors of absorption and the radiation pressure for porous particles, and investigate the influence 
on the y8-ratio. Blowout sizes are calculated for various stellar luminosities and porosities, and an approximation equation is derived 
to estimate the blowout size as a function of these parameters. Furthermore, we investigate the influence of the porosity on the dust 
equilibrium temperature. 

Conclusions. The blowout size increases with the particle porosity and stellar luminosity. The dust temperature of porous particles is 
lower than the one of the compact spheres, in particular the temperature of blowout grains decreases for porous particles. 

Key words, stars: circumstellar matter - protoplanetary disks - infrared: stars ~ interplanetary medium - scattering 



1. Introduction 

Debris disks contain micrometer-sized dust grains, which can be 
observed in scattered light in the optical to near-infrared wave- 
length range and by thermal emission at longer wavelengths. 
To interpret observations of debris disks these particles are of- 
ten assumed to be compact, homogenous, and spherical (e.g 
IRoccatagliata et al.||2009| |Ertel et al.||20lT| |2012j). This is a 
convenient simplification, since it makes the Mie scattering for- 
malism applicable ( |Mie||1908] |Wolf & Voshchinnikovl|2004[ l. 
However, there are various processes that would tend to create 
more complex, e.g. inhomogeneous or almost arbitrarily shaped, 
particles, instead. In cold clouds gas molecules can freeze out 
on the surface of dust grains, so that the particles are com- 
posed of multilayer, inhomogeneous materials ( |Krugel|2003[ ). In 
young, gas-rich circumstellar disks, collisions of particles occur 
in which the dust grains either agglomerate and form larger par- 
ticles or get shattered in fragments with complex shapes. In this 
context, porosity might support the grain growth process in the 
disk ( |Blum & Wurm|2008| l. 

Modeling recent observations of dust in various environ- 
ments also motivates the concept of porous dust grains. For ex- 
ample, |van Breemen et aL] ( |201 l| l used a distribution of hollow 
spheres to fit the 9.7//m silicate absorption profile of the ISM 
based on the model of |Min et aL] (,2005 ). Demyket al.| ( |1999[ ) ap- 
plied some degree of porosity, which was necessary to reproduce 
silicate features in the spectrum of dust around the protostellar 
objects RAFGL7009S and IRAS 19110-H1045. Scattered light 
images of the circumstellar disk suiTounding the T Tauri star IM 
Lupi may indicate the presence of fluffy aggregates (Pinte et al. 
[20D8) . |B]rns"tiel et al.| ( |20T0l ) and |Ricci et aT] ( |20T0||20T2| ) used 
porous spherical grains to specify the dust opacity in a sample 



of young protoplanetary disks. For a general study of the influ- 
ence of porosity on the dust opacity in protoplanetary disks, see 
ISemenov etaL] ( |2003l ). 

Modeling debris disks is an important part of understand- 
ing the evolution and formation of stellar disks, hence the 



evolution of planetary systems. With that in mind, Augereau 
era r] (fT999l[200T] l, |Wyatt & Dent| ( [2002l l, and 'Fitz gerald et al 



(2007]l used porous grains for disk modeling, resulting in suit- 
able approaches to reproduce observational data of HR 4796 A, 
jSPictoris, Fomalhaut, and AU Microscopii, respectively. Acke| 
|et al. ( |2012| l found evidence of fluffy aggregates in the disk of 
Fomalhaut. These are required for the thermal emission charac- 
teristic of small grains and the scattering behavior of large parti- 
cles. A dependence exists between grain porosity and dust tem- 
perature ( Greenberg & Hage|1990 Voshchinnikov et al.|2006 1, 
causing a shift in the infrared emission to longer wavelengths. 

Our study is motivated by the fact that debris disk model- 
ing using compact grains often reveals a significant discrepancy 
between the minimum particle size and the calculated blowout 
size (e.g. Roccatagliata et al. 2009 Ertel et al. 201 1) . The in- 
fluence of porosity on the radiation pressure has already been 
considered in previous studies. |Mukai et aL] ( |1992| l and |Kimiira| 



et al. (1997 1 examined the j8-ratio for fractal aggregates using 
effective medium theory, while|K6hler et al.|(|2007 ) applied dis- 
crete dipole approximation. We investigate porous particles and 
outline various aspects as the absorption efficiency factor, the 
blowout size, and dust temperatures. In Sect. |2]we provide an 
overview of the applied numerical tool for calculating the effi- 
ciency factors and describe our model of porous dust grains. The 
analysis of the absorption cross sections is presented in Sect. [3] 
and the blowout size as a function of selected parameters is dis- 
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cussed in Sect. |4] Finally, we present the equilibrium tempera- 
ture of porous dust particles in Sect.|5] 

2. Calculation methods 

2.1. Porous grain model 

Within the infinite variety of porous particles we consider those 
that are spherical and from which material is removed so that 
voids are formed. The advantage of this model is the small num- 
ber of parameters describing each particle. A debris disk con- 
tains a vast number of grains of different shapes, sizes, and orien- 
tations. The appearance of the disk results from averaged optical 
properties. This motivates for adopting a grain model that re- 
sults from an ensemble of particles with similar geometric prop- 
erties. The spherical grain model with voids satisfies this task, if 
the voids are distributed uniformly over the volume and thus the 
grain orientation is negligible. 

Porosity is the ratio of the vacuum volume to the total vol- 
ume of the grain. 



^ ~ ^vacuum / ^total 

P e [0, 1] . 



1 - Vsolid/V. 



total; 



(1) 



For example, a compact sphere is given for !P = 0. The radius 
a is denoted as the minimum radius of a sphere that encases the 
grain. 

Two grains with equal radius and porosity can have funda- 
mentally different hole sizes, considering the examples of a hol- 
low sphere and a grain that contains many small holes. In prin- 
ciple, holes of diff'erent sizes can be assumed within a grain. 
Voshchinnikov et al. ( 2005] ) adopted a size distribution of the 
holes in their model, in which the probability of the presence of 
a hole is assumed to be inversely proportional to its volume. The 
holes in our grain model are determined both by the number of 



dipoles and the process of the grain generation (Sect. 2.3 i 



2.2. DDSCAT 

We used the program DDSCAT to calculate the scattering and 
absorption cross sections of irregular particles ( Draine & Flatau 
1994 2010[ l. It is based on the theory of discrete dipole approx 



imation (DDA), where a continuous particle is approximated by 
the corresponding spatial distribution of discrete polarizable 
points on a cubic lattice. The great advantage of the DDA is 
the high flexibility of the particle structure, but the applicability 
is limited by the conditions that the lattice constant dj must be 
small compared to the wavelength and to the scales of the in- 
ternal grain structure. The first condition is roughly expressed 
by (Draine & Goodman 1993 Draine & Flatau 1994 and Draine 
[2000l i 



\n\kdi < 1, 



(2) 



where n is the complex refractive index and k - 27r//l the circular 
wavenumber This criterion is satisfied for |n - 1| < 3 ( |Draine &| 
[Flatau 2010 ). With the volume of the material in the particle, 

V = N(dif = 4/37rfl3(l - P), it follows a = (iz^)'^^ t/i = 
'"WT.'^^^ condition^ 



|n| 27r' 



a < 0.09873 



A I N \^'^ 
M\l-p) ■ 



(3) 



For our studies we use astronomical silicate (|Draine | |2003b|cl l 
with a bulk density of p = 3.5 g cm"^ ( Draine 2003a i. 
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Fig.l. The limit of the applicability of the DDSCAT for 
a compact spherical grain of astronomical silicate, repre- 
sented by ~500,000 dipoles. According to condition |3] the 
DDSCAT delivers suitable results for particles with radius a < 

N " 



tests with the 



0.09873j^(Y^) '~ (solid line). In Appendix 
DDSCAT are presented, and we define the crosses as the upper 
limit. 



Before applying the DDSCAT to investigate the behavior of 
porous grains, we verified the optimum parameter settings (num- 
ber of dipoles used, etc.) in the case of compact spherical grains. 
Therefore, we calculate the absorption cross section with the 
DDSCAT and compare the results with absorption cross sections 
i) calculated via Mie theory and ii) using the DDSCAT with a 
number of dipoles decreased by a factor of 2. The results are 
summarized in Appendix [A] 

In Fig. [T] the applicability of the DDSCAT is illustrated for 
a compact spherical grain, represented by ~500,000 dipoles. 
According to condition |3] the DDSCAT delivers suitable results 
for the area below the solid line. While for the solid line the re- 
fractive index of the silicate is taken into account, the crosses 
are derived by the tests with the compact grains (Appendix [A|. 
These crosses mark the used upper limit of applicability. 

2.3. Generating porous grains 

We briefly describe the process of generating porous model 
grains: A^' - 519,832 dipoles are arranged on a cubic lat- 
tice, forming the shape of a compact sphere. A random num- 
ber generator is used to remove single dipoles to create grains 
with a porosity P. The remaining dipole arrangement contains 
N - (I - P) ■ N' dipoles, forming a sphere with voids. Each 
dipole represents a cubic subvolume (dxf of the surrounding 
material, while a removed dipole represents a cubic void of the 
same size. The large number of dipoles and the application of the 
random generator ensure a statistical distribution of the holes. 
Thus, no preferential direction exists and the optical properties 
are roughly independent of the grain orientation. The number of 
dipoles for the particles used in this paper are listed in Table 
[T] The particle cross sections are shown for different porosi- 
ties in Fig. |2] Some dipoles seem to be detached from the rest 
of the complex, but it should be considered that these dipoles 
may be connected to the particle in the direction perpendicu- 
lar to the presented cross section. Nevertheless, the procedure 
of particle generation does not prevent the occurrence of single 
dipoles, which are isolated from the rest of the configuration. 
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'P=o.6 



Fig. 2. Particle cross sections through the middle of the grain for 
porosities from !P = 0. 1 to !P = 0.6. 



Table 1. Number of dipoles for particles with porosity P. 
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0.1 


0.2 


0.3 


0.4 


0.5 


N 


467849 


415866 


363882 


311899 


259916 



For P - O.l the statistical number of isolated dipoles is negligi- 
ble if compared to the number of dipoles of the entire grain, but 
this number rises sharply with porosity. Therefore, only particles 
with porosities up to !P = 0.6 are considered in this paper. 



2.4. p-ratio and blowout size 

We imagine a spherical grain in the orbit of a central star that 
is attracted by stellar gravitation and repelled by the radiation 
force. With decreasing radius, the radiation force becomes in- 
creasingly important. As a consequence, there is a minimum 
particle radius qbo, called blowout size, at which the particle is 
blown out of the system and which can be calculated as a func- 
tion of the grains' optical properties, the stellar luminosity L», 
and stellar mass M« ( [Burns et aL]|1979| |Artymowicz| 1 988| l. By 
modeling the density and grain size distribution of debris disks, 
a minimum grain radius is derived, which is often significantly 
above (a factor of ~2 up to one magnitude in some cases) the 
value calculated under the assumption of compact grains (e.g. 
HD 104860, HD 8907, HD 377, HD 107146, HD 61005, HD 



191089; IRoccatagliata et al. 2009 Ertel et al. 2011). Despite 



model degeneracies, particularly between minimum grain sizes 
and inner disk radii, there have been several attempts to explain 
this discrepancy. A more detailed grain size distribution, derived 
in collisional evolution studies of debris disks and resulting in 



an overdensity of small particles ( Thebault et al. 2003 Krivov 
|et al . 2006 2008 ), can explain a factor of two. Grain segregation 
([Roc catagliata et al. 2009 ) or the occurrence of ice ( [Greenber^ 
& Hage,,1990^ affect the value of the blowout size as well 
( Augereau et al.|[2001[ l. However, another possible explanation 



for this problem can be the nature of the porous particles, which 
requires more investigations, since porosity influences the scat- 



tering and absorption properties of the particles ( 


Voshchinnikov 


|& Mathis|1999||Voshchinnikov et al.|2005||2006 


1. 



2.4.1 . Expression for the y8-ratio 

We now derive an expression for the yS-ratio, the ratio of radiation 
(f^iad) to gravitational force (Fgra), 



y6: 



^rad 
1 gra 



(4) 



Under the assumption of an optically thin disk, the resulting ra- 
diation force is given by 



Rl 



.1 r- 



(i)B,(r,)d^, 



(5) 



(| Kohler & Mann|2 002 2004 ) where the star is considered as a 
black body with effective temperature r», radius /?» and mass 
M«. The quantity 8^(7^) is the Planck function, c is the speed of 
light, and QpriA) is the wavelength-dependent efficiency factor 
of the radiation pressure. 



Gpi- = Qabs + Gsca(l -g), 



(6) 



which is a function of the absorption (Qabs - Cabs/(^fl^)) and 
scattering cross section (2sca = Csca/i^a^)), and the asymme- 
try factor g. For a given grain mass = 4/37rfl-'(l - 'P)p the 
gravitational force on the particle amounts to 



Mtm„ 



(7) 



where y is the gravitational constant. Using the Eqs.|4j|5]and|7] 
delivers 



1 1 1 r 
^ a{l -r)pj 



with 



117.76 



kgs^ Rl 
~m^ W,' 



(8) 



(9) 



For the Sun one thus derives C* = C© = 2.865 ■ 10"" s^m"^. 
The quantity C» is a function of stellar parameters alone, while 
the other factors in front of the integral in Eq.[8]represent the dust 
properties. The integral reflects the interaction between dust and 
star. In Table |2] the factor is presented for main sequence 
stars of spectral type A to K. The stellar parameters are taken 
from|Binney & Merrifield| ( |199^ , as published in |Schmidt-Kaler| 
( [T982^ 

The particle is ejected from the system, if fiia) > 1/2. A 
particle that initially moves on a circular orbit is compelled on an 
elliptic orbit for^S < 0.5. The trajectory is parabolic (hyperbolic) 



for /3 = 0.5 (> 0.5), and the particle leaves the system ( Krivov 
let al.|2006D . 



Table 2. The physical properties of main sequence stars ( Binney 
i& MeiTifield„1998 ) and respective stellar constant C« (see text 
for details). 



Spectral type 


T, [K] M, [Me] R, [Re] 


C, [Co] 


U [Lo] 


AO 


9520 


2.9 


2.4 


1.99 


42.5 


A5 


8200 


2.0 


1.7 


1.45 


11.7 


FO 


7200 


1.6 


1.5 


1.41 


5.4 


F5 


6440 


1.3 


1.3 


1.3 


2.6 


GO 


6030 


1.05 


1.1 


1.15 


1.4 


G5 


5770 


0.92 


0.92 


0.92 


0.84 


KO 


5250 


0.79 


0.85 


0.91 


0.49 


K5 


4350 


0.67 


0.72 


0.77 


0.17 
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We would like to remind the reader that in the literature 
(see [Burns et al.|[T979l |Artymowiczl[T988l |Krugel|[2003| [Hahii 
2010] ! an approximation equation is usually given that neglects 
the wavelength dependence of 2pr- In the Umit of the geometri- 
cal optics, Qpi is equal to 1, resulting in 



(5 = 0.577 



gcm-3/ \fim) \Lo/\Mo 



(10) 



Equation [To] is easy to handle and therefore suitable for a rough 
estimation of the blowout size. However, for more precise stud- 
ies (as required in our investigations) the wavelength dependence 
of Qpr expressed in Eq.[8]must be considered. 

2.4.2. Previous studies 

The influence of porosity on the radiation pressure has already 
been considered in previous studies. Mukai et al. ( 1992 i,|Kimura| 



larger than the ones for compact spheres. Around A ~ I jjm, the 
curves have a wavy pattern, which diminishes with increasing 
porosity. 

Furthermore, we calculated the absorption cross sections of 
the same grain with an eff'ective medium theory (Fig. |5]l, us- 
ing the Bruggeman rule (Bruggeman |T935]l. Besides the main 
similarity of both figures, the main diff'erences are found for 
A < 0.2 i^m and A > 20 fim. The values for the porosity P = 0.6 
are 0.25 to 1.2 times larger than the ones for compact spheres. In 
addition, the wavy pattern around A ~ I //m is more pronounced 
and can be detected also for !P = 0.6. 



4. Blowout size for different parameter 

We determine the blowout size as a function of the porosity and 
the stellar luminosity. 



etaL|(fT997), and|Kohler etaL](|2007|i have examined porous ag- 4. 1 . Blowout size for porous particles in the solar system 



gregates with fractal structures, ballistic particle-cluster aggre- 
gates (BPCA), and ballistic cluster-cluster aggregates (BCCA), 
which are described in terms of the fractal dimension (e.g. 
|Witten & Cates|1986| l. These aggregates are two extreme cases 
of relatively compact and highly porous aggregates, whereas 
we apply a grain model that allows us to increase the porosity 
step by step, permitting an analysis of the physical dust prop- 
erties with changing porosity. Mukai et al. (1992) investigated 
the jS-ratio based on Maxwell Garnett effective medium the- 
ory ( [Maxwell Garnett|1904ll, and[Kirnura et al.|(|1997|l u s ed the 



Bruggeman rule ( Bruggeman| 1935 1. Hage & Greenberg ( 1990| l 
assessed the accuracy of Qabs when determined with Maxwell 
Garnett, while Bazell & Dwek ( 199 ^ report a discrepancy at 
longer wavelengths. Voshchinnikov & Mathis ( 1999| found that 



effective medium theories are considerably worse for calculat- 
ing scattering cross sections and the asymmetry factor g than for 
extinction cross sections. 

In jKohler et al] ( p006| [2007] ) the DDA is applied to sim- 
ulate light scattering by irregular aggregates. They calculated 
the optical properties of larger grains only for one wavelength 
(A - 0.6 /im) and thus neglect the wavelength dependence of 
2pr. In our study the number of dipoles is much larger, which al- 
lows us to calculate the optical properties for a wide wavelength 
range and for grains up to a = 4 fim. We then use the obtained 
data to investigate not only the /3-ratio, but also the blowout 
size for various porosities, stellar temperatures, and number of 
dipoles, and simultaneously calculate the dust temperature of 
porous grains, especially the one of the blowout grains. This 
gives a comprehensive overview of porous grains in debris disks. 

3. Results for the efficiency factor of absorption 

The j6-ratio of a dust grain is a function of the optical properties. 
The simulation software DDSCAT delivers the quantities Qabs. 
Qsci, and g for each wavelength and radius, therefore we inves- 
tigate the absorption cross sections for porous particles with ra- 
dius a - I /im (upper image of Fig. [3]l. The porosity varies in 
steps of 0.1 from 'P = 0.0 to 0.6. To emphasize the dependence 
on porosity, the curves are normalized to the values of compact 
spheres (Qabs.o) in the bottom image of Fig. |3] The relative dif- 
ferences of Qabs increase with porosity for all wavelengths lower 
than 25 fim. Qabs drops with increasing porosity for wavelengths 
from 0.2 to 25 fim and increases for A < 0.2 /im. For /I > 50 /im, 
Qabs increases with porosity up to !P = 0.3 and then decreases 
again. The values for the porosity P = 0.6 are 0.21 to 1.13 times 



At first we investigate the impact of grain porosity on the 
blowout size of the solar system. For this purpose, the /3- 
distributions of compact spheres and porous particles are deter- 
mined. The effective temperature of the Sun is r» - 5777 K, and 
the porosity varies in steps of 0.1 from T = 0.0 to 0.6. The dis- 
tributions are presented in Fig.|4j and the corresponding blowout 
sizes are listed in Table |3] 

1 . The j6- value strongly increases with radius, reaching a maxi- 
mum, and then starts falling off. With increasing porosity, the 
maximum of the yS-distribution decreases and shifts towards 
larger radii. Thus, the maximum is the largest for compact 
spheres. Furthermore, the slope at the larger radii becomes 
flatter 

2. The blowout size of compact spheres is abd = 0.46 //m. For 
comparison, the value 0.33 fim is obtained with the approxi- 
mation Eq. 10 In Mukai et al. ( 1992| ) a blowout size between 

0. 3 fim and 0.35 fim can be estimated, and in Kohler & Mann[ 
(2002) between 0.45 fim and 0.5 fim (see paper. Figs. 4b and 

1, resp.). 

3. The higher the porosity, the larger the blowout size. For P = 
0.6 the blowout size is abo = 0.84//m, 1.8 times larger than 
for T = 0.0. This factor is comparable to the factor between 
the minimum grain size obtained with disk modeling and the 
value calculated under the assumption of compact spheres. 

4.2. Blowout size of stars of various spectral type 

The influence of different stellar luminosities on the blowout size 
is investigated. For this purpose, we consider compact spherical 
grains with particle radii up to 4- fim, since higher luminosities 
increase the radiation force and thus the blowout size. The values 
for Tt and C» for the eight different spectral types from KO to 
A5 are summarized in Table |2] The resulting yS-distributions are 
shown in Fig. [6] 

1 . The higher the luminosity, the higher the yS-values, whereby 
the blowout size increases. In addition, the maximum of the 



Table 3. Blowout size as a function of grain porosity P in the 
solar system. 



■p 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


cbo [j"rn] 


0.455 


0.500 


0.549 


0.603 


0.678 


0.744 
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Fig. 3. Wavelength-dependent absorption cross section Qabs as a 
function of porosity V. The particle radius is a = 1.0 jum. Note 
the two characteristic silicate peaks at /i ~ lOjum and 18/^m. 
In the bottom image the results are normalized to the values of 
compact spheres, Qabs.o- 
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Fig. 5. 
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Same as Fig. [3 
Bruggeman ru 



but derived with effective medium theory 
e). 



0.5 
a [|j.m] 

Fig. 4. The y6-ratio as a function of the grain radius a and porosity 
V for the solar system. The intersections of the different distri- 
butions with the 0.5 limit mark the blowout size abo- 



Temperature 




Fig. 6. The y6-ratio as a function of the grain radius a for stars of 
various spectral type (different stellar temperature). The particles 
are compact spheres. 
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jS-distribution shifts to smaller particle radii for higher lumi- 
nosities. 

2. For some luminosities one finds /? < 0.5 for all particle radii, 
so there no blowout size exists in these systems. In particu- 
lar, only stars with temperatures r» > 5250 K have a blowout 
size. There is also no blowout size at lower temperatures for 
porous grains, since the maximum of the y6-distiibution is the 
greatest for compact spheres, and it decreases with increas- 



ing porosity (Sect. 4. 1 



3. The blowout size at r» = 9520 K is larger than 4/im, 
and even higher values are expected for porous particles. 
Calculations of the optical properties of such large grains 
exceed our computational capabilities, therefore stars with 

> 9520 K are not investigated further. 

4. Several temperature curves intersect the 0.5 limit in the ex- 
amined particle radii interval even twice, whereby two val- 
ues for the blowout size exist (e.g. the = 5770 K curve). 
Further forces are acting on very small particles, such as the 
Poynting-Robertson effect and the Lorentz force, which are 
negligible for larger grains, but have to be taken into account 
for the balance of forces on very small grains (Griin et al. 
[2001) 1. 

4.3. Parameter space 

As shown in the previous two sections, the blowout size depends 
on both the grain porosity and the stellar temperature. We now 
determine the blowout size for each parameter configuration. 
The computational effort is large, especially particles with radii 
a > 2 /urn increase the calculation time dramatically. We gener- 
ate and evaluate the data for the following parameter space: 



0.2, 0.3, 0.4, 0.5, and 0.6, 
= 5770 K, 6030 K, 6440 K, 



- seven porosities: T = 0.0, 0.1 

- five stellar temperatures: r» 
7200 K, and 8200 K, 

- six number of dipoles, from 54, 432 to 5 19, 832 dipoles 



We vary the total number of dipoles and thus the size of the inner 
voids to investigate the influence of different particle structures. 
The calculated blowout sizes of the special particle presented in 
Sect. 2.3 are illustrated in Fig. |7] and listed in Table |4] Missing 
table entries imply a blowout size larger than 4 yum and are thus 
not computable within the parameter settings. 

Figure [T] again reveals the blowout size increasing both with 
grain porosity and stellar temperature. For r» - 5770 K the 
blowout size of porous grains CP - 0.5) is larger by a factor 
of 1.5 if compared to compact spheres. The higher the tempera- 
ture, the stronger the effect, so that the factor amounts to 2.8 for 
Tt - 7200 K. These factors are close to matching the existing 
discrepancy mentioned in Sect. 2.4 The variation in the number 
of dipoles causes a dispersion of the blowout sizes. This disper- 
sion increases also both with !P and T,. In most cases, abo rises 
with decreasing number of dipoles. However, there are excep- 
tions. The largest deviations from the other calculated blowout 
sizes occur for the lowest number of dipoles, the coarsest parti- 
cles, whose flBO are significantly larger than those from the other 
ones (triple dotted lines). 

Figure |7] can be used to estimate the blowout size for a 
given grain porosity !P and stellar temperature Tt. The differ- 
ent blowout sizes at constant porosity and temperature, caused 
by the varying number of dipoles, set the range of applicability 
in which the blowout size can be found. 
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Fig. 7. Blowout size cibo as a function of porosity V and stellar 
temperature for various number of dipoles. Different colors 
represent different temperatures. The values for the special par- 



ticle presented in Sect. 2.3 are illustrated as the solid line. 



Table 4. Blowout size in //m as a function of grain porosity V 
and stellar temperature r» . 















p 


5770 


6030 


6440 


7200 


8200 


0.0 


0.42 


0.60 


0.83 


1.29 


2.15 


0.1 


0.45 


0.67 


0.94 


1.50 


2.45 


0.2 


0.50 


0.75 


1.07 


1.75 


2.94 


0.3 


0.54 


0.86 


1.25 


2.12 


3.67 


0.4 


0.60 


0.99 


1.52 


2.69 




0.5 


0.64 


1.18 


1.93 


3.64 




0.6 


0.70 


1.51 


2.75 







4.4. Fitting model 

In this section an approximation function f{'P,T„) for the 
blowout size is derived. We find that the results of the special 
particle presented in Sect. 2.3 can be analytically approximated 
very well using the following equation: 



flBO i'P, n) = (^1^-1-^2) ■ + ihT, + h) , 



ki = 5.07 ■ 10"^ — , ^2 = -28.85 yum, 
K 

k3 = 7.01 ■ lO""* — , ^4 = -3.65 yum, 
K 



(11) 



a= 1.781- 



61.72 K 
- 5670K' 



It contains seven coefficients (three for a, and the kj, i e N<4). 

The fit is presented in Fig. [8] The relative deviation between 
fit and measured values is for the highest temperature lower than 
8 %, so Eq. 11 is a good way to estimate the size of the blowout. 
An online tool for calculating the blowout size based on this 
approximation function is available at 
http://wwwl.astrophysik.uni-kiel.de/blowout/. 

4.5. /3 influenced by other effects 

We discuss other possible effects which might influence the /3- 
ratio. 
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Fig. 8. Calculated blowout sizes (solid line) and derived approx- 
imation function (dotted lines). Using the equation, it is possible 
to estimate abo as a function of grain porosity and stellar tem- 
perature. 



4.5.3. Ejection limit 

The ejection limit ySejection = 0.5 only applies to particles that 
have been produced by planetesimals that move in circular or- 
bits, i.e. for which the orbit eccentricity is e = 0.0. The general fi 
threshold depends on the position of the elliptical orbit, where 

the particle is set free, and is ySo.ejection = \ ( i + ' cos(O) )- 
centricities as low as e = 0.1 can cause /3-values from 0.45 to 
0.55. Owing to the different ejection conditions for each parti- 
cle, the blowout size depends on the eccentricity. This influence 
is quite strong compared to the other effects discussed above. 



4.5.4. Furtlier forces acting on dust grains in debris disks 

The derivation of Eq. |8] for fi only considers gravitational and 
radiation forces. In addition, the Poynting-Robertson effect and 
stellar wind drags are acting on the particles, and in the case of 
charged particles also the Lorentz force. The magnitude of these 
forces depends strongly on the particle size. The gravitational 
force dominates for particles with radii a > I jjm. For particles 
with radii 0.05 fim or lower, the Lorentz force prevails. For in- 
termediate particle radii, the radiation force is important beside 
gravitation and Lorentz force. In|Grun ( 1999^ , the various forces 
are listed and compared for different particle radii. 



4.5.1. Particle modeling 

The particles used in this paper have a spherical shape, into 
which holes are formed. Besides prolate and oblate shapes, real 
particles can have a much more complex non symmetric struc- 
ture. However, the vast variety of particle shapes cannot be ex- 
amined. The spherical shape with inclusions was chosen as rep- 
resentative of all porous and arbitrary shaped particles. 

It should be noted that large particles are generally more 
porous than smaller ones, because they can be seen as agglomer- 
ates of smaller subunits. Lamy et al. ( 1987j l modeled the density 
of small grains by p = 2.2 - 1.4 (a^ + a), where a is in jjm and p 
in gcm"^. With increasing particle radius this density law leads 
to higher j0-values. 

There is no unique definition of the optical properties and 
bulk densities of what is considered as an "astronomical sili- 
cate". As the y6-ratio is inverse to the density, the choice of p has 
a strong impact on the blowout size. It is also conceivable that 
individual particles consist of several materials. 

The mentioned points do not lead to errors in the calcula- 
tions. Instead, they show the possible complexity and variety 
of particles in the disk, hence the difficulty in specifying the 
blowout size of a stellar system, since a separate y6-value can be 
allocated to each particle and even to each particle orientation. 



4.5.2. Central star 

The stellar component is considered as a black body with the 
effective temperature T,. In contrast to this, Lamy (1974i and 
|Artym owicz 



1988) motivated to use an observed or a modeled 
photospheric spectrum (e.g. Kurucz model). We did not consider 
this in our parameter space study because fundamentals are first 
investigated. In the treatment of a real object, the shape of the 
spectrum should be used. Furthermore, C, is calculated for pa- 
rameters of a main-sequence star model. Real stars can deviate 
from these values, so real values should be used for stellar mass 
and radius if available. 



5. Dust temperature of fluffy grains 

Finally, the influence of porosity on the dust equilibrium tem- 
perature is investigated. For this purpose the dust temperature 
Tg is determined if the grains are only heated by direct stellar ra- 
diation. When taking the energy conservation into account, the 
distance r of the particle from the star is 



2" ^ (rg) gabs (a, A) dA 



(12) 



( Wolf & Hillenbrand|2003[ ). We calculated the dust temperature 
at a distance of r = 50 AU for stellar temperatures Tt = 5770 K 
and Tt = 7200 K and illustrate them in Figs. |9] and [TO] respec- 
tively. 

1 . The temperature of the smallest particles increases with the 
grain radius up to a maximum before it falls off steeply for 
larger radii. In addition, the maximum shifts with increasing 
porosity from 0.28 fim for P = 0.0 to 0.5 fim for f - 0.6 in 
the case of = 5770 K, and from 0.24 pm for T = 0.0 to 
0.43 pm for P = 0.6 in the case of = 7200 K. 

2. For particles with radii 0.03 pm < a < 1.2 pm the following 
is valid. The higher the porosity, the colder is the particle. 
For a fixed radius, the drop in temperature decreases with 
increasing porosity. For P = 0.6 and radii a < 0.03 pm, the 
temperature is even higher than for corresponding radii at 
r = 0.5. 

3. The blowout sizes as calculated in Sect. l4.3l are marked in the 
figures and reveal that most of the presented dust particles 
will be blown out of the system. 



The temperature of the blowout grains is an important pa- 
rameter, because this is the quantity that is derived from observa- 
tions. We calculated the dust temperature of the blowout grains 
Tgo at different distances from the star, since it might depend 
on the cooling efficiency, and hence on the stellar distance (Fig. 
[TT) . The temperature of the blowout grains is at all distances for 
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Fig. 9. Dust temperature Tg at the distance of 50 AU from a star 
with the effective temperature r» = 5770 K, depending on the 
grain radius a and porosity !P. The blowout sizes are marked as 
crosses. 
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Fig. 10. Same as Fig . |9j only for the stellar temperature r» 
7200 K. 
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Fig. 11. Dust temperature of the blowout size grains Tbo as a 
function of porosity P and the distance r from the central star 
The stellar temperature is - 5770 K. 
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Fig. 12. Dust temperature Tg at the distance of 10^ stellar radii 
from a star with effective temperature r» = 2500 K, depending 
on the porosity and the effective radius Oeff = (1 - 'P)'^^ a- The 
dust material is a 1:1 mixture of astronomical silicate and amor- 
phous carbon. 



porous grains slightly lower than for compact spheres. For ex- 
ample, at r = 5 AU and for P - 0.6, the dust temperature is 22 K 
bel ow the one for P = 0.0 . At 100 AU this difference is 9.2 K. 

jVoshchinnikov et al. ( |2006) determined the dust temper- 
ature for porous particles using the effective medium theory 
once and the method of layered spheres once. The material was 
a mixture of astronomical silicate ( Laor & Draine|[T993| l and 
amorphous carbon ( |Rouleau & Martin[ 1991 1 with a weight- 
ing of 1:1. Comparing our grain model with the studies from 
IVoshchinnikov et aL] ( |2006| l, we furthermore calculated the dust 
temperature for a mixture of astronomical silicate and amor- 
phous carbon, which is presented in Fig. 12 as a function of the 
effective radius, flgif = (1 - 'P)'^' a. The central star has an effec- 
tive temperature of = 2500 K, and the assumed grain distance 
is r = 10^ R,. 



The results of IVoshchinnikov et al. 



(2006 1 are comparable 
to the ones calculated with the DDA. For the layered spheres 
the temperature decreases with increasing porosity. This state- 
ment applies to the EMT only up to a porosity of T - 0.7, 
when the temperature increases again. The decrease in tem- 
perature for lower porosities can also be observed in our ap- 
proach. Greenberg & Hage (19901 ascertained in their stud- 
ies the increasing of temperature at very high porosities, where 
the temperature approaches the value of subunits of the porous 
particle. Certainly the behavior differs for the lower porosities, 
since the temperature remains constant or increases slightly. 
,Voshchinnikov et aL] ( [2006 j have also performed investigations 
for different stellar temperatures and distances and found the 
same behavior For high porosities, the differences between the 
method of layered spheres and the EMT could reach ~ 15 %, 
and for compact grains they were less than 1 %. 

In conclusion, porous particles with radii up to 1.3yum are 
generally colder than compact spheres, and the maximum of the 
infrared emission is thus shifted to longer wavelengths. Our re- 
sults agree very well with those of Voshchinnikov et al. (2006j). 



6. Conclusions 

When modeling the density and grain size distribution in debris 
disks, the minimum particle size is often significantly larger than 
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the corresponding blowout size. We investigated porous dust 
grains, which provide a possible explanation. Since the intro- 
duced y6-ratio depends on the optical properties of the grain, we 
determined the absorption cross sections of porous particles as a 
function of the porosity. The blowout sizes of six different sam- 
ple stars and seven porosities (including compact spheres) were 
calculated. We found that the blowout size increases with grain 
porosity and stellar temperature. It should be noted that both the 
choice of the bulk density, as well as the particle ejection limit, 
has a strong influence on the obtained results. An approxima- 
tion equation was derived from the numerically calculated val- 
ues with which the blowout size can be estimated. Finally, the 
dependence of the dust temperature on particle radius and poros- 
ity was investigated. It was found that porous particles are colder 
than compact spheres, and the temperature of blowout grains in 
particular decreases for porous particles. 

The optical properties of porous particles were determined 
with the DDSCAT for an extended range of wavelengths and 
particle radii. These efficiency factors of absorption, scattering, 
and radiation pressure are summarized in several tables and are 
available for further applications. An online tool for calculating 
the blowout size based on the derived approximation function is 
available at http://wwwl.astrophysik.uni-kiel.de/blowout/. 
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Appendix A: Tests with compact grains 

A.1. Comparison Mie theory and DDA 



We use the program mi ex ( |Wolf & Voshchinnikov|2004 1 to cal- 
culate the absorption cross sections based on JVlie scattering and 
compare them with DDA results. Eight grain radii from 0.5 to 
lOOjum and 200 wavelengths from 0.1 to 1000//m are exam- 
ined (Fig. 13 I. With an increasing size parameter, the DDSCAT 
computing time increases, and the results become highly in- 
accurate, depending on the chosen number of dipoles, in our 
case 519,832. In Fig. A.l the relative deviations of the re- 
sults of DDA (2abs) and Mie (Qabs, Mie) are presented, A^abs = 
l(2abs - Gabs,]Mie)/2abs,]Miel- One Can scc in Fig.Othat at short 
wavelengths the DDA values are significantly below the Mie 
ones. Thus these results are neglected, exceeding the limits of 
the DDA, and there is a lower limit of the wavelengths A^i„ for 
which at all wavelengths A > /Imi,, follows: Aq^b^iA) < 5%. Here, 
/l,ni„ increases with the grain size a (Table A.l i. The limit of ap- 



plicability can be raised by increasing the number of dipoles, 
resulting in an increase in the calculation time as well. 



0.25 




1 10 100 1000 1 10 100 1000 

X [p.m] X [nm] 



Fig. 13. Absorption cross section Qabs derived with DDA 
(crosses) and Mie (solid lines) as a function of the wavelength. 
The particles are compact spheres with radii a from 0.5 to 5.0 fim 
(above) and 10.0 to 100.0//m (bottom). 



Fig. A.l. Relative deviations of the absorption cross sections de- 
rived with DDA and Mie, A^abs = l(2abs - 2abs,Mie)/2abs,Mie I- 

The particles are compact spheres with radii a from 0.5 to 5.0 jum 
(above) and 10.0 to lOO.Oyum (bottom). 
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In conclusion, to calculate absorption coefficients for wave- 
lengtfis sliorter than 0.5 fim, the DDA can be applied within our 
parameter settings for grains with radii a < 5 fim. Despite these 
limitations, the similarities of both approaches are remarkably 
good. 



Table A.l. The smallest wavelength /Imin, at which the relative 
deviations of the DDA and Mie results are not larger than 5 %. 
The number of dipoles is 519,832. 



a [fim] 

/IminljUm] 


0.5 1.0 2.0 5.0 
<0.1 0.11 0.20 0.44 


a [/um] 

'ImmljUm] 


10.0 20.0 50.0 100.0 
0.87 1.7 4.2 ~ 12 



A.2. Decreased number of dipoles 

The convergence of DDA depends on the particle shape, and 
for compact spheres it is easier to reach convergence (|Yurkin 
|et al.|[20T0l l. To check the accuracy of the method we decrease 
the number of dipoles by a factor of 2 and compare the results. 
The decreased number of dipoles causes only small changes. The 
lower limit of the wavelengths /Imin is slightly increased (Table 
A.2 1. These values mark the upper limit of applicability of the 
DDSCAT (Fig.[T). Besides that, the chan ges are negligible. 



Table A.2. Same as Table 



A.l 



but for the relative deviations 



between normal and halved number of dipoles. 



a [//m] 
/l„,iiil;^m] 


0.5 1.0 2.0 5.0 
<0.1 0.12 0.24 0.55 


a \jim] 
/Imml/'m] 


10.0 20.0 50.0 100.0 
1.08 2.2 5.1 ~ 13 



References 

Acke, B., Min, M., Dominik, C, et al. 2012, Astronomy and Astrophysics, 540, 
A125 

Artymowicz, P. 1988, The Astrophysical Journal, 335, L79 

Augereau, J. C, Lagrange, A. M., Mouillet, D., Papaloizou, J. C. B., & Grorod, 

P. A. 1999, Astronomy and Astrophysics, 348, 557 
Augereau, J. C, Nelson, R. P., Lagrange, A. M., Papaloizou, J. C. B., & Mouillet, 

D. 2001, Astronomy and Astrophysics, 370, 447 
Bazell, D. & Dwek, E. 1990, The Astrophysical Journal, 360, 142 
Binney, J. & Merrifield, M. 1998, Galactic Astronomy, Princeton, NJ : Princeton 

University Press, (Princeton series in astrophysics) 
Bimstiel, T, Ricci, L., Trotta, F., et al. 2010, Astronomy and Astrophysics, 516, 

L14 

Blum, J. & Wurm, G. 2008, Annual Review of Astronomy & Astrophysics, 46, 
21 

Bruggeman, D. A. G. 1935, Annalen der Physik, 416, 636 
Burns, J. A., Lamy, R L., & Soter, S. 1979, Icarus, 40, 1 

Demyk, K., Jones, A. R, Dartois, E., Cox, R, & D'Hendecourt, L. 1999, 

Astronomy and Astrophysics, 349, 267 
Draine, B. T. 2000, The Discrete Dipole Approximation for Light Scattering by 

Irregular Targets, ed. Mishchenko, M. I., Hovenier, J. W., & Travis, L. D., 131 
Draine, B. T. 2003a, Annual Review of Astronomy and Astrophysics, 41, 241 
Draine, B. T. 2003b, The Astrophysical Journal, 598, 1017 
Draine, B. T. 2003c, The Astrophysical Journal, 598, 1026 
Draine, B. T. & Flatau, P. J. 1994, Journal of the Optical Society of America A, 

11, 1491 

Draine, B. T. & Flatau, R J. 2010, ArXiv e-prints 

Draine, B. T. & Goodman, J. 1993, The Astrophysical Journal, 405, 685 
Ertel, S., Wolf, S., Metchev, S., et al. 2011, Astronomy and Astrophysics, 533, 
A132 



Ertel, S., Wolf, S., & Rodmann, J. 2012, Astronomy and Astrophysics, 544, A61 
Fitzgerald, M. R, Kalas, R G., Duchene, G., Pinte, C., & Graham, J. R. 2007, 

The Astrophysical Journal, 670, 536 
Greenberg, J. M. & Hage. J. I. 1990, The Astrophysical Journal, 361, 260 
Griin, E. 1999, Interplanetary Dust and the Zodiacal Cloud, 673-696 
Griin, E., Gustafson, B., Dermott, S., & Fechtig, H. E. 2001. Interplanetary Dust, 

Astronomy and Astrophysics Library (Springer), 804 
Hage, J. I. & Greenberg, J. M. 1990, The Astrophysical Journal, 361, 251 
Hahn. J. M. 2010, The Astrophysical Journal, 719. 1699 

Kimura, H., Ishimoto, H., & Mukai, T. 1997, Astronomy and Astrophysics, 326, 
263 

Kohler, M., Kimura, H., & Mann, I. 2006, Astronomy and Astrophysics, 448, 
395 

Kohler, M. & Mann, I. 2002, Proceedings of Asteroids, Comets, Meteors - ACM 
2002., 500, 771 

Kohler, M. & Mann, I. 2004, Journal of Quantitative Spectroscopy & Radiative 
Transfer, 89, 453 

Kohler, M., Minato, T., Kimura, H., & Mann, I. 2007, Advances in Space 
Research, 40, 266 

Krivov, A. V., Lohne, T, & Sremcevic, M. 2006, Astronomy and Astrophysics, 
455, 509 

Krivov, A. V., Miiller, S., Lohne, T. & Mutschke, H. 2008, The Astrophysical 
Journal, 687, 608 

Kriigel, E. 2003, The physics of interstellar dust, loP Series in astronomy and 

astrophysics, Bristol, UK: The Institute of Physics 
Lamy, P. L. 1974, Astronomy & Astropysics, 35, 197 

Lamy, P. L., Gruen, E., & Perrin, J. M. 1987, Astronomy & Astropysics, 187, 
767 

Laor, A. & Draine, B. T. 1993, The Astrophysical Journal, 402, 441 

Maxwell Garnett, J. C. 1904, Royal Society of London Philosophical 

Transactions Series A, 203, 385 
Mie, G. 1908, Annalen der Physik, 330. 377 

Min, M., Hovenier, J. W, & de Koter, A. 2005, Astronomy and Astrophysics, 
432, 909 

Mukai, T, Ishimoto, H., Kozasa, T, Blum, J., & Greenberg, J. M. 1992, Annual 

Review of Astronomy & Astrophysics, 262, 315 
Pinte, C, Padgett, D. L., Menard, F., et al. 2008, Astronomy and Astrophysics, 

489, 633 

Ricci, L., Testi, L., Natta, A., et al. 2010, Astronomy and Astrophysics, 512, A15 
Ricci, L., Trotta, F., Testi, L., et al. 2012, Astronomy and Astrophysics, 540, A6 
Roccatagliata, V., Henning, T, Wolf, S., et al. 2009, Astronomy and 

Astrophysics, 497, 409 
Rouleau, F & Martin, P G. 1991, The Astrophysical Journal, 377, 526 
Schmidt-Kaler, T. 1982, Landolt-Bornstein Numerical Data and Functional 

Relationships and Technology, 2 
Semenov, D., Henning, T, Helling, C, Ilgner, M., & Sedlmayr, E. 2003, 

Astronomy and Astrophysics, 410, 611 
Thebault, P., Augereau, J. C, & Beust, H. 2003, Astronomy and Astrophysics, 

408, 775 

van Breemen, J. M., Min, M., Chiar, J. E., et al. 2011, Astronomy and 

Astrophysics, 526, A152 
Voshchinnikov, N. V., Il'in, V. B., & Henning, T. 2005, Astronomy and 

Astrophysics, 429 

Voshchinnikov, N. V., Il'in, V. B., Henning, T, & Dubkova, D. N. 2006, 

Astronomy and Astrophysics, 445, 167 
Voshchinnikov. N. V. & Mathis, J. S. 1999, The Astrophysical Journal, 526, 257 
Witten, T. A. & Gates, M. E. 1986, Science, 232, 1607 
Wolf S. & Hillenbrand, L. A. 2003, The Astrophysical Journal, 596, 603 
Wolf, S. & Voshchinnikov, N. V. 2004, Computer Physics Communications, 162, 

113 

Wyatt, M. C. & Dent, W. R. F. 2002, Monthly Notices of the Royal Astronomical 
Society, 334, 589 

Yurkin, M. A., Min. M., & Hoekstra, A. G. 2010. Physical Review E, 82, 036703 



10 



